in 

o ; P. Strack 

<N 

>v 
a ■ 

Kinetics of Oscillating Neutrinos 

> 

in 
o 
in 
o 
in 
o 

pu,: 

> 

A Master's Thesis Submitted to the Faculty of the 
DEPARTMENT OF PHYSICS 
In Partial Fulfillment of the Requirements For the Degree of 
MASTER OF SCIENCE 
In the Graduate College 
The University of Arizona 
(Accepted and defended April 2005) 



Copyright © Philipp Strack 2005 



Statement of the Author 



This Master's thesis has been submitted in partial fulfillment of requirements for an advanced de- 
gree at The University of Arizona and is deposited in the University Library to be made available 
to borrowers under the rules of the Library. 

Brief quotations from this Master's thesis are allowable without special permission, provided 
that an accurate acknowledgment of source is made. Requests for permission for extended quota- 
tion from or reproduction of this manuscript in whole or in part may be granted by the copyright 
holder. 



SIGNED: Philipp Struck 



Acknowledgements 



I thank the selection committee of the Akademische Auslandsamt der Universitdt Heidelberg 
for providing me with this wonderful opportunity to study at the Departments of Astronomy 
and Physics of the University of Arizona. The Heidelberg professors Dr. J. Stachel and Dr. C. 
Wetterich have helped me with the application for my studies in the USA and I thank them. I 
acknowledge financial and organizational support from the German Fulbright-Commission, and 
the Baden-Wuerttemberg Stiftung. They have made life easier. 

I am indebted to my supervisor Prof. Dr. Adam Burrows for intense 1 and superb scientific 
guidance. This collaboration has been an honor for me. 

The main part of this thesis, the generalized Boltzmann formalism in fiat spacetime is pub- 
lished in Physical Review D (Strack & Burrows 2005). We acknowledge support for this work 
through the SciDAC program of the Department of Energy under grant number DE-FC02- 
01ER41 184. In particular, we thank Prof. Dr. Z. Chacko for reading the manuscript of the paper 
and for useful discussions, and C. Meakin and M. Prescher for helpful conversations and for read- 
ing the manuscript of this thesis. The latter has also provided invaluable help with technicalities 
and latex and has thus softened the author's inability to use computers. 

The members of my thesis committee consisting of the professors Dr. A. Burrows, Dr. Z. 
Chacko and Dr. J. P. Rutherfoord, gave generously of their time to help me for which I thank 
them wholeheartedly. I express my gratitude to Prof. Dr. W. J. Duschl who has supported me 
more than once with applications. 

This work was carried out at Steward Observatory and I gratefully acknowledge the kind 
hospitality I have received. The graduate students here at Steward Observatory have been very 



1 Efforts leading to this thesis were started at the beginning of September 2004 and finalized by the end of 
March 2005 and were at all times accompanied by three graduate courses in physics. 



6 Acknowledgements 



integrative and as representatives for all of them I wish to thank R. S. Bussmann, I. Momcheva, 
H. J. Seo, N. Siegler, S. Sivanandam, A. Diamond-Stanic, and J. Trump. 

My whole family and in particular my father, Rolf Dieter Strack, and my beloved partner in 
life, Anne Gerrit Knepel, have showered me with love. I hope I can return this one day. Gerrit 
has also proof-read the whole thesis and I thank her. 



Dedicated to my grandfather Georg Walter Strack (1923 -) 



Abstract 



In the standard approach to neutrino transport in the simulation of core-collapse supernovae, one 
calculates the spatial, temporal, and spectral evolution for each neutrino species in terms of the 
classical density in phase-space or the associated specific intensity. The neutrino radiation is cou- 
pled to matter by source and sink terms on the "right-hand-side" of the transport equation and 
together with the equations of hydrodynamics this set of coupled partial differential equations 
describes, in principle, the evolution of core collapse and explosion. However, with the possibil- 
ity of neutrino oscillations between species, a purely quantum-physical effect, how to generalize 
this set of classical Boltzmann equations to reflect oscillation physics has not been clear. To date, 
the formalisms developed have retained the character of quantum operator physics and have 
not been suitable for easy incorporation into standard supernova codes. In this thesis, we de- 
rive generalized Boltzmann equations for quasi-classical, real-valued phase-space densities with 
which any numerical code for the classical transport equation can easily be generalized to include 
neutrino oscillations. This new, extended set of equations incorporates the full mixing Hamilto- 
nian consisting of the vacuum oscillations, the matter-induced effective mass for the v e , and the 
self-interactions from neutrino-neutrino scattering. The equations are intrinsically nonlinear for 
two reasons: blocking correction at high densities and the off-diagonal v-v Hamiltonian which 
maintains SU(N) flavor symmetry for an N-flavor system. Standard oscillation phenomenology 
including resonant flavor conversion (the MSW effect) and the interplay between decohering 
matter coupling and flavor oscillations are retained in a quantum-physically consistent fashion. 
The extension of the formalism to curved spacetime is discussed. One has to take into account 
two modifications. Firstly, spinors in curved spacetime must be analyzed with the spin connec- 
tion in a noncoordinate basis and, secondly, for the kinetic terms on the left-hand side of the 
Boltzmann equation the Christoffel connection must be employed. Both changes are brought to 
an algebraic footing via vielbein (tetrad) fields associated with the square root of the metric. The 
generalized kinetic equations in curved space derived in this work describe neutrino evolution in 
the most general setting encountered in astrophysical contexts such as the early universe or in 
the vicinity of supernova cores. 
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Introduction 



Mixed particles and their oscillations are fundamental for a wide range of interesting physics: 
quark-mixing by the Cabbibo-Kobayashi-Maskawa matrix, its leptonic analog for massive neu- 
trinos (Bilenky & Petcov 1977 and 1987), hypothetical photon-axion and photon-graviton os- 
cillations in the presence of external magnetic fields (Moertsell, Bergstroem & Goobar 2002, 
Raffelt & Stodolsky 1988) and K° - K° oscillations (Akhmedov, Barroso & Keraenen 2001). 

Physically, these quantum systems are coupled to macroscopic systems and through external 
interaction their quantum evolution is altered. One prominent astrophysical context in which such 
oscillation for macroscopic systems is important, involves the neutrinos in supernova cores that 
may execute flavor oscillations while simultaneously interacting with ambient supernova matter 
(Takahashi, Sato, Burrows & Thompson 2003, Thompson, Burrows & Horvath 2000, Burrows 
& Gandhi 1995, Raffelt 1996, Pantaleone 1995, Qian & Fuller 1995, and Wolfenstein 1979). 

The primary motivation of this work is to provide a straightforward generalization of the 
Boltzmann formalism with which to analyze the kinetics of oscillating neutrinos with collisions 
including energy redistribution and inelastic scattering. By taking ensemble-averaged matrix el- 
ements of quantum field operators for mixed particles, following the pioneering work in the 
references (Raffelt & Sigl 1993, Raffelt, Sigl & Stodolsky 1992, Akhiezer & Peletminskii 1981, 
and de Groot, van Leeuwen & van Weert 1980), we obtain quasi-classical phase-space densi- 
ties that satisfy simple Boltzmann equations with coupling terms that account for the neutrino 
oscillations. The formalism is clear, numerically tractable, and does not contain operators or 
wavefunctions. 

Finally, we discuss how to extend our formalism to curved spacetime. This includes the 
treatment of generalized resonances similar to the MSW effect. We present kinetic equations for 
oscillating neutrinos interacting with a background medium in curved spacetime. 
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The plan of this thesis is as follows. 

In chapter 12 we introduce the Wigner phase-space density operator approach from which 
we derive our formalism involving classical phase-space neutrino flavor densities and their 
off-diagonal, overlap correlates. The latter couple the different flavor states to account for 
neutrino oscillations. 

In chapter[3] the different contributions to the mixing Hamiltonian originating in the vacuum 
oscillations, matter-induced effective mass for the v e , and the neutrino-neutrino scattering 
are discussed. The full set of Boltzmann equations is derived for both the flavor phase-space 
densities and the associated specific intensities. The advantage of the specific intensity rep- 
resentation is that collision terms can be written in a simple and intuitive fashion. 

In chapter^ we demonstrate the usefulness and applicability of our formalism by success- 
fully testing it on its quantum-mechanical accuracy. We demonstrate with several simple 
examples that the set of equations reproduces what one would expect from the standard quan- 
tum mechanics wavefunction approach. 

Firstly, we give analytic solutions to the standard vacuum flavor oscillations for neutrinos 
in a box. Secondly, we consider flavor oscillations with absorptive matter coupling. This is 
frequently referred to as quantum decoherence (Raffelt, Sigl & Stodolsky 1993). Thirdly, 
it is shown that our formalism is capable of calculating the neutrino evolution in phase- 
space when the resonant matter-induced flavor conversion (MSW effect) is important. We 
show complete analytic and numerical congruence to the existing wavefunction formalism as 
discussed in the papers by Wolfenstein 1978, Bethe 1986, Mikheyev & Smirnov 1986 . 
We show that our formalism can provide a tool to analyze active-sterile neutrino conversion 
in a beam with absorptive matter coupling. One can observe "indirect decoherence" of the 
active-sterile conversion. This means that even though the sterile neutrinos do not interact 
with a medium their oscillatory pattern is indirectly influenced since the v e 's execute absorp- 
tive collisions. 

In part|nj we sketch how to make the relevant changes if one wishes to apply our formalism 
in a curved spacetime. This is potentially important in astrophysical contexts such as the early 
universe and in the vicinity of supernova cores. 

In chapter [6j we analyze in a heuristic fashion how neutrino oscillation phenomenology is 
modified in curved space. Therefore, we introduce the vielbein formalism which is the nec- 
essary ingredient to do calculations involving spinors in the framework of the general theory 
of relativity. Furthermore the vielbein fields are useful tools to write the Liouville operator 
on the left-hand-side of the Boltzmann equation in a manifestly covariant fashion to account 
for the curvature. This is done in section lTTl 
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• Culminating in our efforts, we present new kinetic equations for oscillating neutrinos with 
collisions in curved spacetime in chapter^ The equations of section l731 include the possi- 
bility of spin-flips into a sterile right-handed species by coupling of the anomalous magnetic 
moment of the neutrinos to an external magnetic field. In principle, eqs. (I7.24> govern the 
evolution of the quasi-classical density in phase-space for interacting, massive neutrinos in- 
cluding all known effects - from the quantum-mechanical viewpoint as well as from the 
gravitational perspective. 



Parti 



Flat Spacetime 
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Wigner Phase-Space Density Operator 



2.1 Construction 

To analyze a multi-particle system and reflect its inherent statistics, one usually sets up the den- 
sity matrix of the system (Akhiezer & Peletminskii 1981, de Groot, van Leeuwen & van Weert 
1980). For classical systems this is the phase-space density. For quantum fields, the conceptually 
most similar analog is the Wigner phase-space density operator (Wigner 1932 & 1984). In what 
follows we shall briefly sketch how to construct this quantity. 

For Dirac neutrinos, we begin with the Dirac equation 

i(A-m)|J = 0, (2.1) 
which describes the dynamics of spin-^ fields, iff, in Minkowski space-time with the metric 

rf = diag (1,-1, -1,-1) . (2.2) 
y are Dirac matrices, which satisfy 

{y,7 v } = 2/T- (23) 

Curly brackets denote matrix anti-commutators. Following the standard second quantization pro- 
cedure, we write for the momentum eigenstate expansion of general Dirac fermion annihilation 
operator in position space: 

<Kr, f) = J ijL £ a^ip, E p )e-* r e- iE >> + 2>*V(p, E v )e* r e iE *' , (2.4) 
and for the Dirac fermion creation operator in position space: 
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^(r, t) = j J] *;V(p, £ p ) e -'' pr e -' £ "' + apV(p, E^j**' , (2.5) 

where u s (p, t) and v ? (p, f) are four-component Dirac spinors describing fermionic and an- 
tifermionic states, respectively. The summation is over additional quantum numbers, e.g. spin 
states of the particle. We now tailor if/ir, t) to the physical properties of the neutrino and define 

<A(r, t)=f |^ [flp(0"(p. E P ) + il p (f)v(p, £ p )] e' pr . (2.6) 



A Fourier transformation connects momentum and position space representation of the oper- 
ators: 

fl (p,r)= f -^3M t (p,£ P M(r,0e-* r . (2.7) 

Note that here the operators are explicitly time-dependent. This corresponds to the Heisenberg 
picture. For mixed neutrinos, the representation in Fourier space cannot diagonalize the energy 
dependence and the momentum dependence simultaneously. 

The Dirac spinors «(p, E p ) and v(p, E p ) refer to negative-helicity particles and positive- 
helicity antiparticles, respectively. 

We need to incorporate the appropriate spin-statistics properties into the description. At equal 
times the Pauli exclusion principle is implemented by anti-commutation relations of the creation 
and annihilation operators in momentum space: 

Mp, f), a V, 0) = MP, t), & V, f)) = (2n?S 3 (p - p') . (2.8) 

One can now write the Wigner phase-space density operator for Dirac neutrinos in terms of the 
creation and annihilation operators in momentum space as 

p(r, p, t) = j l^e-'P'Vtp - Ip', t)a(p + Ip', t) . (2.9) 

For antineutrinos we have the expression 

p(r, p, t) = J |^e-*'V(p - Ip', t)b(p + Ip', . (2.10) 

The same Wigner phase-space density operator can be formulated in terms of the creation and 
annihilation operators in position space as long as the integration variable is changed: 

P(r, p,f) = J 7^3 e ~ , ' Pr > t ( r - 5 r '> *Mr + Ir*. - (2-11) 
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and similarly for antineutrinos 



P(r, p, t) = j -^L-e-^'fir - l -r' , t)fa + \r' , t) . (2. 12) 



The Wigner phase-space density operator possesses the sought-after properties of a quantum- 
mechanical generalization of the classical phase-space density. To see this, we integrate eq. (12. 1 li 
with respect to p and use the representation of the delta-function in Fourier space to write 

j </ 3 pp(r,P,0 = l<A(r,f)l 2 , (2.13) 

which resembles the correct observable probability. When eq. ( 12.1 1> is integrated with respect to 
r one has 

j d 3 rp(r,p,t) = 

which is again the correct expression. This can be achieved by a shift of variables: 

v 

u , (2.15) 

and a subsequent relabeling of the integration variable. The above derivation was done for 
fermions. The procedure is identical for bosons however. Then, we impose commutation rela- 
tions instead of the anti-commutation relations for the operators. Furthermore no restrictions to 
the mass of the mixed particle of interest are made. 

In the general case of mixed particles (this includes oscillating fermions such as neutrinos but 
also mixing bosons such as photon-axion conversion systems), the creation and annihilation op- 
erators naturally yield mixed particles. To project these composite objects onto their"irreducible" 
subspaces we take the matrix elements in the appropriate basis. We simultaneously perform an 
ensemble average. We show that this transforms the Wigner phase-space operator into its matrix 
elements. The matrix elements are of great use since the diagonal elements are quasi-classical 
phase-space densities. In the off-diagonal of this matrix we have macroscopic overlap functions. 
One is now calculating numbers and has not to worry about the time evolution of operators over 
Fock space. Instead the value of a number changes with time. 

In principle, these matrix elements are very similar to the density matrix or statistical operator 
of a quantum-mechanical system as for instance for the spin density matrix of particle beams. 



c/ 3 nA(r, t)e~ ipr 



I<A(P,0I 2 



(2.14) 
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The quantum number spin can formally be brought on the same page as the quantum number 
flavor (Friedland & Lunardini 2003). The description of both systems respects the SU(N) sym- 
metry for a N-composite problem. The form of the interaction is alike and has been useful as a 
calculational analogy. 

In the following section 12.21 we demonstrate how the Wigner phase-space operator yields 
quasi-classical phase-space densities. 
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2.2 Matrix elements 

To decompose the Wigner phase-space density operator into quasi-classical phase-space densi- 
ties, scalar functions with seven-dimensional domain, we need to calculate its matrix elements. 
For clarity, we set the dimensions of flavor space to two. We average over an ensemble consisting 
of two mixing particle species and take the matrix elements in the number-density basis of Fock 
space. This results in the expression 

r={ni\p\nj) = ^:^y (2 - 16) 

where the indices i, j run over e, fi, and * means complex conjugation 1 . For the off-diagonal 
macroscopic overlap functions we write 

fen = <«vXr,p,0|nv„> 

= J §^ eip ' r <«v.l flt (p - jp'. + jp'- f > k) ' < 2 - 2 °) 

and for the v e neutrino phase-space density we have the expression: 
f Ve = (n Yt \p(r,-p,t)\n Yt ) 

-^e ip ' r (n Ve | fl f (p - ip', t)a(p + ip', r) \n Ve ) , (2.21) 
t) j 1 2 2 1 

and identically for the v'^s. 




1 Following (de Groot, van Leeuwen & van Weert 1980), one can define the covariant energy-momentum 
densities 

r""(x) = c J ^jfp^T (2.17) 
and the conserved neutrino current densities 

-Lfxx\r-f \ , (2.18) 

wherein, by conservation of matter, 

3^=0. (2.19) 

x? = (ct,x) denotes a position four- vector and p* = (E/c, p) the four- momentum respectively. With the 
aid of these entities, eq. 12. 181 and eq. 12. 171 which describe a quasi-classical field, one can apply the 
framework of classical field theory to investigate for instance the symmetries and conserved quantities 
of the neutrino ensemble. 
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Heuristic ally, we can put this process of simultaneous projection and averaging as follows. 
The integration process is symmetric with respect to p'. Thus, each of the a'(p-ip', f)<z(p+5P', 
hits every momentum value of the n Vr . Regarded this way, we can anti-commute a'(p - jp', f) 
and a(p + 4p', f) in momentum space until they are adjacent and possess identical arguments. At 
that instant they act on their eigenstates and yield numbers as their eigenvalues. Subsequently, we 
are left with real-valued functions possessing seven-dimensional domain; quasi-classical phase- 
space densities f Vj . 

For completeness, we list the matrix elements for the associated antiparticles: 



f = (hi\p\hj) = 




(2.22) 



Note that it is the hermiticity of p and p that implies the off-diagonal terms of its matrix ele- 
ments be complex-conjugates of one another. 

A generalization to more particle species is straightforward. The diagonal terms are real- 
valued and denote quasi-classical phase-space densities. The off-diagonal entries are complex- 
valued macroscopic overlap functions. They represent the degree of coherence or macroscopic 
overlap of the ensemble under consideration. For completely decohered ensembles and for non- 
mixing ensembles the off-diagonal entries vanish. More precisely, f e ^(r, p, f) for instance gives 
the number of v/s that result from incoherently adding up the quantum mechanical amplitudes 
from each otherwise v M . This requires that flavor-entangled many -body states can be neglected 
or are small in number. 



2.3 Conclusions 

• The salient results of chapter[5]are eqs. (I2.16> and (I2.22> . 

• We have abandoned the operator character of the mathematical entities describing the neu- 
trinos, eq. i2.9\ . and are now in position to work with real-valued phase-space densities and 
their complex-number-valued overlap correlates, eq. (12. 161 . 

• In the following, we construct dynamical equations for the quasi-classical phase-space den- 
sities which are very similar to the equations of classical transport theory (see for instance 
eq. (13341 ). 
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Generalized Boltzmann Equations 



3.1 2x2 Mixing Hamiltonian 

Naively, we can derive a dynamical equation for the matrix elements of the Wigner phase-space 
density operator by merging the Boltzmann equation with the Heisenberg equation. A derivation 
from first-principle quantum mechanics has been performed in many different contexts (Landau 
1957, Akhiezer & Peletminskii 1981, Rudszky 1990, Raffelt, Sigl & Stodolsky 1993). We need 
to account for the matrix character of the phase-space densities and, following the nice papers 
of Raffelt & Sigl 1993, generalize the "Liouville terms" on the left-hand side by introducing 
anti-commutators resulting in the expression 

or i r dr\ i f. ar\ . rri _ „ 

(3.1) 

Here, the collision matrix is given by 

C(r,p,f) = L Cv ' C r v ") , (3.2) 



where the off-diagonal elements describe flavor-changing collisions. In the standard model, flavor 
number is conserved at every vertex. So by restricting to these interactions one can in general set 
the off-diagonal elements to zero. 

Oscillation phenomenology is incorporated via the commutator on the right-hand-side, the 
term that originated in the Heisenberg equation, with the mixing Hamiltonian given by, 

Q(r, p) = Q v . dC (s) + Onat(r) + &w(r, p, t) - Q„ v {v, p, t) . (3.3) 

In this cumulative expression, £3 mc (s) is the usual vacuum contribution where mixing comes 
about through the non-degeneracy of masses and the entailed distinction between flavor or weak 
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interaction eigenstates and the mass eigenstates. £? mat (r) is the ordinary matter contribution that 
most notably results the effective mass for the v e neutrino in matter. Q vv is the non-diagonal 
Hamiltonian that describes the neutrino self-interactions or neutrino-neutrino scattering and Q V y 
is the non-diagonal neutrino-antineutrino interaction Hamiltonian. 

Making use of the antineutrino phase-space densities as given in eq. (I2.22i the dynamical 
equation for the anti-particles looks completely analogous, 

— — - + - { v, — — } + - <p, — — } - -i\Q,T\+C , 
dt 2 \ <9r J 2 \ dp J L ' J 

(3.4) 

where apart from the obvious change in the collision terms and the different phase-space densities 
(denoted by a tilde), we must reverse sign for the coupling coefficients in Q mdt , Q vv and Qy V , 

Q(r, p) = Q mc (s) - Qmat(r) - i2„(p) + Qy V (p) . (3.5) 

The vacuum mixing contribution is alike for neutrinos and its antiparticle, 

tic I — cos 2.0 sin 2$\ 
^vac(e)=— • r, , (3.6) 



L \ sin 26 cos 20 

where 9 vac is the mixing angle for vacuum neutrino oscillations between v e 's and v^'s and L is 
the vacuum neutrino oscillation length: 

Antics 

L = —5-4 , (3.7) 

where s is the neutrino energy, mi and mi are the masses of the neutrino mass eigenstates, and 
Am 1 -m\ - m\. The other variables have their standard meanings. 

The electron neutrino finds the most leptonic interaction partners in the ambient matter. 
Therefore, we have the matrix: 

Ona<(r) = ( 2 q °J , (3.8) 

where A is the dimensionless neutrino-matter interaction amplitude (Mohapatra & Pal 2003, 
Giunti 2004): 

2 V2G F 

A(r) = — — F -n e (r) , (3.9) 



where Gf denotes Fermi's constant, and n e the electron number density. The interaction ampli- 
tude A, through n e , is a function of spatial position. For antineutrinos the sign of A is reversed as 
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Fig. 3.1. Charged-current scattering reaction leading to the effective mass of the v e in an electron back- 
ground. W is the charged current gauge boson mediating the weak force. 



indicated in eq. (13.51 . 

Diagrammatically, the charged-current reaction contributing to the effective mass of the v e is 
given in Fig. d3.1> . 

Clearly, the v^'s and v r 's find orders of magnitude less scattering partners in the ambient 
matter. Arguments can be given that the neutral-current interactions involving the Z boson are 
the same for all flavors and thus are not included in the mixing Hamiltonian; here we are only 
interested in differences between the flavor species (Mohapatra & Pal 2003). For an illustrative 
scenario where this neutral-current scattering can be important see section f7"2l and Fig. 17. II 

3.2 Neutrino self-interactions 

The neutrino evolution in a neutrino background is nontrivial and continues to raise considerable 
discussions in the literature (Pantaleone 1992 A, Pantaleone 1995, Qian & Fuller 1995, Friedland 
& Lunardini 2003, Bell, Rawlinson & Sawyer 2003, Boyanovsky & Ho 2004, Sawyer 2004, and 
Sawyer 2005). Here we sketch the considerations that led to Q vv and work with the most estab- 
lished version that can be found in the literature. 

New is, that we match the wavefunctions to the quasi-classical phase-space densities em- 
ployed in our ensemble-averaged approach. In this work, we do not engage in efforts about 
a possible speed-up of flavor equilibration as proposed by Sawyer, et al.. We retain the one- 
particle description neglecting any effects arising from coherent many-body state formation or 
flavor entanglement (Friedland & Lunardini 2003). The exact expression for the neutral-current 
neutrino-neutrino interaction Hamiltonian for a test neutrino with initial momentum p and the 
background neutrino with initial momentum q is given by (Raffelt & Sigl 1993): 
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H 



1 

2 \ 2 cos 6w 



J dp' J dq' (2tt) 3 <J (3) (p + q - p' - q') ^/^D^y^ 



(3.10) 



where the ^'s are the Dirac neutrino spinors as specified in eq. ( 12. 6L y are the standard Dirac 
matrices as given in eq. J2.3> . cos 6w is the Weinberg angle, g is the weak coupling constant given 
by 



(4 V2G F ) 



1/2 



where Gf is Fermi's constant. D^ v is the Z-boson propagator 



D%(P) 



'z ' 



2 2 

— p 



(3.11) 



(3.12) 



Force mediating gauge bosons that propagate from one neutrino at a given spacetime point to 
another neutrino at another spacetime point is a non-local effect. This can be readily discarded 
since in supernova cores and in most other situations found in nature one is typically interested 
in neutrinos with energies of the order of MeV. This is well below the mass or energy scale of 
the Z-boson of order 100 GeV. Therefore, we neglect the non-local gauge boson effects and work 
with the low-energy local four-fermion coupling. This effective Hamiltonian must satisfy SU(N) 
flavor symmetry for a system of N flavors (Pantaleone 1992 A) and we write the following terms: 



Q vv — 



V2 



v J 



(3.13) 



where the sum is over all neutrino flavors. In this form the SU(N) flavor symmetry is manifest. 
For a two-flavor system consisting of v e 's and v^'s one can rewrite this Hamiltonian to accentuate 
its off-diagonal character (Pantaleone 1992 B, Friedland & Lunardini 2003): 



=fi(l -costfi" 1 ) 



\k\ 2 



(3.14) 



where the coupling coefficient prefers diametral scattering by inclusion of the angle 6 between 
the test neutrino with momentum p and the background neutrino with momentum q. 
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Ve V Vf, / 

y e r ji 

Fig. 3.2. Diagonal contribution to the neutrino-selfinteractions: v e - v e and - scattering. 




Fig. 3.3. Off-diagonal contribution to the neutrino-self-interactions. The Z-boson carries the difference in 
momenta. The graph on the right-hand-side is the low-energy local four-fermion interaction. 



The coupling strength (3 and the normalization of the background neutrino fields are deter- 
mined by the formulae: 



J dV^lf + |<f = 1 



P=^. (3.15) 
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In Fig. 13.21 the two diagonal contributions to the v-v interaction are shown. The time axis 
is vertical. In Fig. 13.31 the off-diagonal contributions to the v-v interaction are depicted. Here 
the point is that in forward direction the momenta of the two scattering partners are exchanged. 
Since the momentum eigenstates are the propagation eigenstates in the ambient matter, this is 
considered as a flavor-exchanging reaction. 

We now match the the v-v Hamiltonian from the wavefunction formalism to our ensemble- 
averaged description. To this end, we make use of our quasi-classical phase-space densities and 
the macroscopic overlap functions of eq. J2. 16b . This involves integration over all the momenta 
q of the ensemble (note that the neutrino ensemble also constitutes the background). 

The v-v mixing Hamiltonian for test neutrinos with momentum p in ensemble-averaged form 
is denoted by 



where, while throwing away the overall phase term in eq. ( I3.14t . which is proportional to the 

identity matrix, we write for the diagonal elements: 



where the momentum integration goes over all momenta in the ensemble. For the off-diagonal 
elements we write in our notation, in terms of the macroscopic overlap functions, 




(3.16) 





(3.17) 




(3.18) 



For the antineutrinos in the ensemble and the ambient matter, 




(3.19) 



And in complete analogy, we write for the interactions of the test neutrino with momentum p 
with the antineutrinos, 
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B eft (p, r,t) = B 



BU 9 ,r,f)=B 



B v Jp,r,t) = B 



B Ve (p,r,t)=B 




d 3 q(l-cos^)/ Ve (q,r,0 
d 3 q(l-cos0 p,1 )/ yf( (q,r,O 
J 3 q(l -cos6» pq )/^(q,r,f) 
A(l-cos0 p «)il(q,r,O 



(3.20) 



The coupling coefficient in front of the integrals has to be implemented in the equations for neu- 
trinos with reversed sign. For a similar way of parameterizing the neutrino-neutrino Hamiltonian 
see the nice paper by Fuller & Qian 1995. 

3.3 Kinetic equations 

We are interested in an analysis of the macroscopic, long-range transport behavior of eq. (13 . II 
and execute the following simplifying procedure. 

• We do not take into account the effects of small-scale external homogeneities on the mixing 
in the ensemble (Sawyer 1990). 

• We average out any small-scale correlations and do not include any other off-diagonal con- 
tributions from the complete Hamiltonian. 

• We neglect off-diagonal contributions from the left-hand-side of eq. ( 13. 11 1. 

• We componentize eq. (13 - 1 1 and include all the contributions from the mixing Hamiltonian in 



To this end it turns out to be useful to define the real part of the off-diagonal macroscopic 



eq. d3~3l . 



overlap in eq. ( I2.16> as 




(3.21) 



and the corresponding imaginary part as 




(3.22) 



respectively. For antineutrinos one proceeds analogously. For the real part one has 
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fr = \(U+r el ) , (3.23) 
and for the corresponding imaginary part there is 

ft ~ 2j {feii ~ fefij ■ 

(3.24) 

With this choice of parameters, everything becomes real-valued and the i's are eliminated 
from the equations. In the standard weak interaction basis or flavor basis the collision terms 
are diagonal and we will not consider flavor-changing reactions here. Henceforth we apply the 
diagonal collision matrix, 

C = f C -°). (3.25) 



. C ,, 

For antineutrinos the collision terms differ and are denoted by 

(3.26, 

Having made these definitions, we are in a position to componentize eq. J3.ll > and eq. 13.41 . 
We work with the general mixing Hamiltonian eq. J3.3i and eq. J3.5I . For two mixing neutrino 
species, here denoted by v e and Vp, interacting with a background medium the generalized Boltz- 
mann equations in their most generic form are: 
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df Ve 9fy e . 8fy e 

— - + v • — - + p • — - 

dt dr F <9p 



-fi | — sin 26 + 2/3 



df*. df. 



df Vu 



dt dr dp 



J (l-cos^)(/ r +/ ( -)j 3 qJ + 
2/3 % j (1 - cos flW) (/; + /;•) rf 3 q + C Ve 
f. sin 20 + 2/3 J (1 - cos flW) ( /r + f r ) d 3 qj - 
2/3 % J (1 - cos (fi + /;•) rf 3 q + C v „ 

i +v -i +p -i^[- (A - cos20)+ 

/3 J (1 - cos?") (/ Ve - / v . - / V( , + 4) rf 3 q] + 
(/v. - / J /» J (1 - cos p(1 ) (/I - rf 3 q 

f +V 'f =(^-^)[f Sin26+/3 f (l-cosPi)(/ r -/^ 3 q - 

/ r [^ (A -cos 20) + 

J*(l - cosPi) (f Ve - / Ve - 4 + 4) J 3 q] . (3.27) 

In a very similar fashion we can write the corresponding antiparticle analog. We need to inter- 
change and add tildes, sign-reverse the matter interaction amplitude A, and complete our set of 
equations: 
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dfi e dfi e . dfi^ 
dt +V ' dr + P dp 



dfv, 54 . df Vfl line . 



-fi sin 20 + 2/3 J (1 - cos 0P«) (/, + /,) d 3 qj - 
2/3 fi J (1 - cos pq ) (J + fi) d\ + C Ve 

126 + 2/3 J (1 - cos 0*) (fi + fi) d 3 qj - 



2/3 fi cos fl«) (/• + ./;■) d\ + C v , 



<3/> 3/r . df r *\2nc , „ s 

i +v 'i +p -^^[- ( - A - cos2e)+ 



dt V ' <9r P <9p 



J*(l - cosPi) (/ y , - / Ve - 4 + 4) d 3 q] + 

(k - k)p j - cos 0Pq ) (/« - /') rf3 q 

: (/v. - /vj ™ sin 26 +fi J (1 - cosei) (/, - f r ) d\ 



~v2nc 

f r [— (-A-cos26) + 



p J (1 - cose') (/ y , - /„, - 4 + / V J d 3 q] . (3.28) 

Neutrino and antineutrino evolution are nonlinearly coupled. The collision terms differ for neutri- 
nos and antineutrinos. Our new generalized eight Boltzmann equations are entirely real-valued, 
yet they contain all the quantum-mechanical oscillation phenomenology. 

This is the complete set of kinetic equations for the real-valued neutrino phase-space densities 
f Ve and f v and the corresponding antineutrino phase-space densities f Vr and f v 1 . In principle, one 
needs to solve the above set of eight real-valued equations to simulate and understand neutrino 
evolution in phase-space as for instance in numerical calculations for the supernova explosion 
mechanism and the early universe. Our generalized Boltzmann equations are intrinsically non- 
linear due to the neutrino self-interactions from Q vv and Q VY . 

Note that in the absence of collisions, Liouville's theorem for the total neutrino flavor specific 
intensity, 

|(/v f +4) = 0, (3.29) 
where 4 denotes the total time derivative, is recovered. 

at ' 



1 For a slightly different way of writing the generalized Boltzmann equations (in terms of the B's in eq. 
13. 16l ~). we refer to the accompanying paper (Strack & Burrows 2005). 
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Moreover, at sufficiently high densities, one must include blocking corrections for the colli- 
sion terms as well as for the oscillation terms which would further increase the nonlinearity of 
the equations. We feel that this is straightforward; there is no need to include these statistical 
properties in the formal framework. 

For blocking corrections to be important for neutrinos, however, the densities must be suffi- 
ciently high (p > 10 8 g cm -3 ). At these densities observable neutrino oscillations are suppressed 
due to the very small effective mixing angle in matter which can be written as (Wolfenstein 1979, 
Pantaleone 1989): 



sin 26 = 



sin 26,,, 



^sin 2 26 vac + (cos26 vac -^ I ) 2 



(3.30) 



In this expression for the effective mixing angle in matter, 6 vac is the mixing angle for vacuum 
neutrino oscillations and in numbers we have for the matter-induced effective mass term: 

As _ 2 yf2G F (Y e /m„)ps 
Am 2 



Am 2 



1.53 x 10 



-2 



gem 



MeV 



10~ 5 eV 2 



Am 2 



(3.31) 



where p denotes the mass density, Y e the number of electrons per nucleon, and m„ the nucleon 
mass. 

One must include the other constituents of the ambient matter fields such as neutrinos, nu- 
cleons and other leptons for an exact effective mixing angle (Qian & Fuller 1995). The density 
dependence is alike and if anything, this additional contribution amplifies the effect of matter 
suppression. 



3.3.1 Transport equations for the neutrino radiation field 

It is instructive to rewrite eq. J3.27b in terms of the specific intensities of the neutrino radiation 
field. This way, the collision integrals can be cast in a simple form. We neglect the neutrino self- 
interactions in the following, having shown how to include them in eq. (13.271 . Also we focus 
on the neutrino degrees of freedom since without the self-interactions neutrino and antineutrino 
evolution are coupled only in the normal fashion through pair processes. We use the one-to-one 
relation between the invariant phase-space densities, f Ve , f v and the specific intensities, I Ve , I v 
(Burrows, et al. 2000, Mihalas 1999) and define 
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The off-diagonal macroscopic overlap terms have the same units as the phase-space densities and 
thus we may define their associated specific intensities as 



.1 en 



'■" (Inhfr 
(2nH) 3 c z 



These physically observable variables are commonly used in hydrodynamic supernova simu- 
lations and we can write the extended set of transport equations including neutrino oscillation 
phenomenology in terms of the specific intensities. 

The generalized Boltzmann equations in the laboratory (Eulerian frame) for the radiation 
field of two oscillating neutrino species with collisions are then: 



2tt 



ldl v v dl v s 3 p d(l Ve e 3 ) 
c dt c dr c dp L 



Idly, v S/ v s 3 p d(l v ^- 3 ) 



c dt 


i 

C 


dr 


- i 

C 


dp 




V 




e p 

c 


d{K^) 


7T- + 

C dt 
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dr 


dp 


ldl ell 


V 


dlen 


e 3 p 
+ — — • 

c 




TT- + 

c dt 


c 


dr 


dp 



2n 

— left sin 29 + C Vjj 
2n 

-— (cos26»-A)J fjU 



2n(I Ve -Iv u \ 
— - sin 29 + (cos 20 - A) 



L 



(3.34) 



These equations and eqs. J3.27b are our major results. The collision terms can be conveniently 
written as (Burrows, et al. 2000) 

C' = -< A + < [ Bv ' + ~r~ f n')I v Xn')dQ' , (3.35) 



1 - n. 



where <t> Ve is a phase function for scattering into the beam integrated over the solid angle dQ' . 
Furthermore, k? v is the sum of all absorption processes n/cr", where n, is the number density 
of matter species i and <x" denotes the absorption cross sections (for scattering processes the 
superscript a is replaced with s). 

is the equilibrium Fermi-Dirac occupation probability 

exp[(e Vf - (jU e - yt/))ySJ + 1 

where fj. e and fi denote the electron chemical potential and the difference between neutron and 
proton chemical potential, respectively, and B Vc is the corresponding black body specific inten- 
sity: 



3.4 Conclusions 
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By. = 




<p-eq 



(3.37) 



Changing the subscript from v e to yields the corresponding parameters for the v^'s. For sterile 
neutrinos, one substitutes v s for and sets the scattering and absorption terms to zero. Similarly, 
one can write a set of equations for antineutrinos with different collision terms and with the 
sign of A reversed. Here, neutrino and antineutrino evolution are implicitly coupled through pair 
processes. 

3.4 Conclusions 

• We have formulated the most generic form of the generalized Boltzmann equations including 
the nonlinear evolution in a neutrino background in eq. J3.27I . 

• We have succeeded to incorporate the purely quantum-physical phenomenon of flavor oscil- 
lations into the classical transport equation (13.341 . 

• In the following, we employ the simple version of our new kinetic equations in terms of the 
specific intensities, eq. (I3.34> . and perform explicit analytic and numerical tests of this for- 
malism in chapter^] 



We show that our formalism is completely consistent with the existing single-particle wave- 
function approach in subsection l4.3.1l 



Simple Tests of the New Formalism 



In practice, it is interesting how the neutrino ensemble behaves with time when there are flavor 
oscillations and collisions simultaneously. Since the interaction rates for the differently flavored 
neutrinos can differ dramatically-depending on the surrounding matter constituents-, the possi- 
bility of flavor oscillations must be taken into account for. 

Eqs. ( I3.27l > and d3.34t serve this purpose. For the illustrative examples in the present chap- 
ter we neglect the neutrino self-interactions, work with eq. ( I3.34L and demonstrate the correct 
limiting behavior of the formalism. 

4.1 Flavor oscillations with absorptive matter coupling 

As a first example, we solve eqs. J3.34i for v e - v M oscillations in a box of isotropic neutri- 
nos that can also experience decohering absorption on matter. This is realized by including the 
neutrino absorption reaction on nucleons as given in Appendix [X] Every absorptive collision 
decoheres the oscillation cycle. In reality, however, due to the small cross-sections of neutrinos, 
frequent collisions with the ambient matter fields require high densities. As indicated above in 
eq. (I3.30> . observable flavor oscillations are suppressed at the high densities found in stellar col- 
lapse (Wolfenstein 1979). To remedy this problem artificially, we "turn off" matter suppression 
by setting the matter interaction amplitude A of Q ma , in eqs. (I3.8i and (13.301 to zero. This is not 
the situation found in nature. 

The goal of the following definitions is to simplify eq. J3.34i for the present purpose and 
write it in terms of dimensionless variables. Therefore, we define the approximate oscillation 
time, 
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and artificially set the interaction rate of the v e 's equal to the interaction rate of the v^'s And 
we define the characteristic absorption time, 



1 (l-K) 

4*c cN A pY„ai e „ 



^ COl "* »T /7 > (4-2) 



where Na denotes Avogadro's number, Y n is the neutron fraction per nucleon, and <x" „ is the 
cross section for absorption on neutrons as given in AppendixlAl The absorption time is flavor- 
independent in our toy calculation. As a "critical parameter" for the characteristic time evolution 
of the system we define the ratio of the oscillation to the absorption timescales: 



a = "vT . (4-3) 

col 

and finally we define the dimensionless time coordinate 

r = -i- . (4.4) 

*osc 

Furthermore, we set B Ve = B Vfi = B v and denote the dimensionless specific intensities by 

%,(r) - \ . (4.5) 
and for the off-diagonal macroscopic overlap functions we write 

a v 
I, 



I e ,(r) = . (4.6) 



Note that the black-body function is flavor-independent, since we take both neutrino flavors to 
be of the same energy (chemical potentials have also been set to zero; they do not affect the 
solutions significantly). The resulting dimensionless version of eq. ( I3.34l > reads: 



1 This definition of flavor-independent interaction rates is again not so realistic. Cross sections depend on 
the mass of the associated leptons which differ by a factor of ss 200. Furthermore the abundance of the 
potential interaction partners in the ambient matter is normally much greater for the v e 's. 
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3/v, 
dr 

m, 



-Jen sin 28 + a (l - 7 Vf ) 
I eM sin 28 + a (l - 7 V(l J 



--T ej[i cos 20 



<5J CU Iy. ~ I 



dr 



- sin 28 + % ejl cos 28 . 

(4.7) 



Since here we consider isotropic neutrinos, it is easy to show that the scattering in and out of 
the beam - represented by the integral over the phase-function - and the scattering part of the 
extinction coefficient cancel 2 . 

In matter and for the initial conditions 7 y = fi ell - I eft = and I Ve + 0, one can derive an 
harmonic oscillator equation for the early rate of evolution of I Vr : 

d% /I 2 



+ - - oT \Iy = const . (4.9) 



3t 2 \2 

As expected, for u« 1, the early time dependence of the solution is predominantly sinusoidal: 
neutrino oscillations dominate. For a » 1, an exponential decay/increase dominates with a 
timescale of I /a. Further algebraic manipulations convert eq. ( 14. 11 into an inhomogeneous ordi- 
nary differential equation of fourth order for the dimensionless electron-neutrino specific inten- 
sity, valid for all t: 



2 In a vacuum (a = 0) and with the initial conditions Z Ve = 1, I v = (and consequently no off-diagonal 
overlap terms at t = 0), we obtain 



X t (j) = 


2 i T 

1 - sin" 28 sin" - 
2 






Mt) = 


sin 2 28 sin 2 - 
2 






K,(T) = 


- sin 29 cos 29 (cos r 
2 




1 


Z,(T) = 


- sin 29 sin t . 
2 



This behavior of the radiation field is unambiguously identical to the probability density obtained by 
squaring the amplitude of a single-neutrino wavefunction in a beam. The off-diagonal terms representing 
the macroscopic overlap peak when mixing of v e and Vf, neutrinos is maximal and vanish when the 
ensemble is single-flavored. 
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d 4 Iy, 



+ 2a 




(a 2 + cot 2 20+l) 




+ a(2cot 2 29+l)—^ + 
v 'or 

(ffcot26») 2 (7^ - l) = . 



(4.10) 



In the following, we evaluate eqs. J4.7i and J4. lQi for three different regimes (without matter 
suppression since the matter interaction amplitude for the v e , A, has been set to zero artificially). 
Initial conditions are I Y = = I ejl = and I Ve = 0.8 and from the large-mixing-angle solution 
(LMA) (Fukuda, et al. "2002) we apply sin 20 = 0.9 and Am 2 = 6.9 x 10~ 5 eV 2 . 

4.1.1 Oscillation-dominated regime: a « 1 

We tailor eq. J4. 101 according to the smallness of a and write 



where we have dropped all terms with a. This is an approximation - even in this regime. 

The time dependence of the solution to this equation is predominantly sinusoidal with oscil- 
lation periods of the order of unity. 

As depicted in Fig. 14. II flavor oscillations dominate the temporal evolution of the neutrinos. 
Asymptotically, the diagonal specific intensities for oscillating v e and neutrinos equilibrate at 
1. 

Absorption on neutrons, and by detailed balance, the resulting emissivity, drive the system 
to radiative equilibrium at the blackbody intensity. The real part of the off-diagonal term, 9i efl , 
takes negative values, whereas the imaginary part, I ejl , oscillates symmetrically around zero. The 
oscillation amplitude of these terms is comparable to the change in magnitude due to absorption. 
Both, Rep and 1 ' efi converge toward zero as the ensemble becomes increasingly decohered; the 
system eventually approaches flavor equilibrium and no oscillations persist. 

In other words, the density matrix of the system becomes diagonal with time due to absorptive 
coupling with matter. 




(4.11) 



4. 1 Flavor oscillations with absorptive matter coupling 
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(a) 
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Time t 



(b) 




10 20 30 40 50 

Time t 



Fig. 4.1. Oscillation-dominated regime, (a): dimensionless specific intensities, (b): macroscopic overlap 
functions. Parameters: e Ve = e Vfl = 10 MeV , p = 4 x 10" g cm 4 , and T = 1 MeV. 
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4.1.2 Oscillating and absorptive regime: a as 1 

In this case, all terms of eq. J4. 10b have to be retained. The solution to the homogenous equation 
consists of a sum of exponentials with the four roots of the following equation in Fourier space 
as exponents: 



where we have set a — I. This equation has two roots of real value and two roots of complex 
value with magnitudes of the order of unity. Thus, decaying exponentials and oscillating expo- 
nentials are represented equally. 

From Fig. 14.21 it is clear that the equilibration due to absorption happens on the same 
timescale as oscillations. The ensemble is guided to flavor and radiative equilibrium. Coherent 
flavor oscillations are disrupted by absorptive collisions. Asymptotically, the diagonal specific 
intensities for the v c 's and v^'s equilibrate. 

Absorption on neutrons, and by detailed balance, the resulting emissivity, drive the v e 's to 
the blackbody intensity. The oscillation amplitude decreases with time; the quantum evolution 
of the system is decohered through absorptive coupling with matter. This is sometimes referred 
to as quantum decoherence (Raffelt, Sigl & Stodolsky 1993). The real part of the off-diagonal 
overlap, ( R ejl , takes predominantly negative values whereas the imaginary part, I efl , oscillates 
symmetrically around zero. Both vanish asymptotically and no oscillations persist. 




(4.12) 



4. 1 Flavor oscillations with absorptive matter coupling 



45 



(a) 




2 4 6 8 10 12 14 

Time t 




Fig. 4.2. v e - Vf, oscillations of isotropic neutrinos in a box with absorptive matter coupling, (a) Specific 
intensities, (b) Off-diagonal macroscopic overlap functions. Parameters: s Ve = e V/J = 10 MeV, p = 8 x 
10 12 gcm- 3 ,T = 5MeV. 
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4.1.3 Absorption-dominated regime: a » 1 

We extract the terms of eq. J4. 10l > consisting of a and write the approximate equation for this 
regime as 

3 3 7v ad 2 T v 18% a cot 2 26- a cot 2 26 

^ + 2^ + 2^ + ^^ 2— =°- (4 ' 13) 

The solution is predominantly exponentially waxing or waning (depending on the initial condi- 
tions) with timescale I /a. 

The equilibration due to absorption is the dominant feature of this regime as shown in Fig. 
14.31 Nonetheless, the specific intensities execute small-amplitude oscillations around the black- 
body value (a small effect due to neutrino mixing). The diagonal specific intensities for v e and 
neutrinos equilibrate at the blackbody intensity in fractions of one reduced oscillation time r. The 
amplitude of the corresponding off-diagonal macroscopic overlap functions is very small (note 
the values on the ordinate). In a nutshell, oscillations do not influence the solution significantly. 



4. 1 Flavor oscillations with absorptive matter coupling 
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Fig. 4.3. Absorption-dominated regime, (a): dimensionless specific intensities, (b): off-diagonal macro- 
scopic overlap functions. Note the small values on the ordinate of (b). Parameters: s Ve = s V(1 = 10 MeV , 
p = 4 x 10 13 g cnT 3 , and T = 10 MeV. 
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4.2 Active-sterile oscillation in a beam with absorptive matter coupling 

The mass of hypothetical sterile neutrinos can only be constrained by observation (Abazajian, 
Fuller & Patel 2001). It is commonly thought that it could be much larger than the mass of 
the active neutrino species resulting in a comparatively large mass-squared-difference: Am 2 ~ 
<9(eV 2 ) (Abazajian, Bell, Fuller & Wong 2004). This shifts the relevance of flavor oscillations to 
relatively high densities. 

As in the preceding section, we artificially "turn off" matter suppression by setting the matter 
interaction amplitude A to zero. We determine the specific intensities of a v e neutrino beam 
simultaneously oscillating into a sterile species v s and being absorbed by nucleons including 
only the reaction in AppendixlAl 

To this end, we define the approximate "reduced" oscillation length 

L 2Hs 
2n sin 26 Am 1 ^ sin 20 



and the characteristic absorption length 



1 (l-gg) 

C N A pY„a^ en ' 



where Na denotes Avogadro's number, and Y n is the neutron fraction per nucleon, and cr^ n is the 
cross section for absorption on neutrons. We define the ratio of oscillation to absorption length 
scales 

(4-16) 

col 



and the dimensionless distance coordinate 



f = r-- (4 - 17) 

lose 



We denote the dimensionless specific intensities and the off-diagonal macroscopic overlap func- 
tions normalized to the blackbody intensity with a hat. 



4.2 Active-sterile oscillation in a beam with absorptive matter coupling 
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(a) 
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Fig. 4.4. Active-sterile conversion with absorptive matter coupling, (a) Specific intensities, (b) Off-diagonal 
macroscopic overlap functions. Parameters: s Ve = s Vs = 40 MeV, p = 5 x 10 13 g crrr 3 , sin 20 = 0.9, T = 10 
MeV, and Am 2 = 0. 1 eV 2 (Abazajian, Bell, Fuller & Wong 2004). 
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The resulting dimensionless version of eq. j3.34l > reads: 



dly, 

dI Vs 



I, 



-I.. + a 





-I es cot 29 



K - h, 

2 



- + < R es co\2Q . 



(4.18) 



Note that there are no collision terms for sterile neutrinos. 

In Fig. 14.41 we depict the solutions to eq. ( 14. 71 for a v e neutrino beam simultaneously oscil- 
lating into a sterile species v s and being absorbed by nucleons. 

In this toy example, a = 0.2. Initial conditions are 7 Vi = fl es = I es = and I Ve — 0.5 . 
Flavor oscillations and collisions happen on the same length scale. While passing through matter 
the beam is guided to flavor equilibrium by the radiative equilibration. The coherence of flavor 
oscillations is decreased with time. Similar to the neutrinos in a box, the v c 's asymptotically 
equilibrate at the blackbody intensity. The sterile neutrinos in the beam converge toward the 
value 0.475 . 

The interacting v e 's indirectly govern the oscillatory pattern of the v,'s even though the sterile 
neutrinos do not interact. The oscillation amplitude decreases with time; the quantum evolution 
of the system is damped through absorptive coupling with matter. The real part of the off-diagonal 
overlap, 7? ej , converges toward the finite negative value -0.35 . The imaginary part, I es , oscillates 
symmetrically around zero and asymptotically vanishes. 



4.3 Matter-enhanced resonant flavor conversion 
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4.3 Matter-enhanced resonant flavor conversion 

When a mixed neutrino traverses matter at certain densities, its flavor composition may be reso- 
nantly altered. This is due to the, compared to the vacuum, large difference between the effective 
neutrino masses in matter. The v e receives an additional effective mass contribution due its (com- 
pared with the other flavors) super-allowed interaction rates. This has been incorporated into the 
equations in this thesis; see eq. J3.8I . 



To demonstrate that our formalism contains the MSW effect (Wolfenstein 1979, Mikheyev & 
Smirnov 1986, Bethe 1986), we solve eq. J3.34i for a monoenergetic one-dimensional neutrino 
beam propagating down a density profile for which resonant matter-enhanced flavor conversion 
takes place. 

We define the dimensionless distance coordinate in terms of the oscillation length 

_ 2tt 

x = x— , (4.19) 

and the dimensionless matter-induced mass term in terms of its resonance value: 

— A A 

A = = . (4.20) 

A res cos 26 

The beam passes the resonance density for A = 1. We analyze the following dimensionless 
version of eq. (I3.34> : 



- I epl sin 26 



I efl sin 29 



dI Ve 
ffx 

ffx 

ffx 

ffx 2 



-cos 20 (l - A) J e/J 

'- sin 26» + cos 26 (l -A)K„ > 

(4.21) 



where C' and C', have been set to zero. 



In Fig. 14. 51 we depict the solutions to eq. ( 14.2 It for a density profile of A = |j . We set ft = 900 
and thus ensure that the scale of spatial inhomogeneities is large compared to the microscopic 
length scales such as the neutrino de Broglie wavelength and the oscillation length. Initially, the 
beam contains only v e neutrinos. The mixing angle is arbitrarily taken to be sin 2 26 - 0.18 (at 
present the LMA is favored (Fukuda, et al. 2002). From Fig. 14. 51 it is clear that at the resonance 
density, ~x = -\^6 = 30, the flavor composition of the beam is radically altered. For higher values 
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Fig. 4.5. MSW effect, (a) Specific intensities, (b) Off-diagonal macroscopic overlap functions. An initially 
y c -beam propagates down the density profile A = Jy, where /? = 900. 
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of x, the beam executes vacuum oscillations. In this illustrative problem the density at production 
is much greater than the resonance density. Then, the spatially averaged survival probability of a 
v e neutrino going from matter to free space should be (Parke 1986) 

(P(v e -> v e )) ~ (1 - P,) sin 2 + P x cos 2 6 , (4.22) 

where P x is the Landau-Zener probability for non-adiabatic transitions. For the chosen density 
profile, propagation is adiabatic and P x = 0. The averaged v e survival probability in Fig. 14.51 
converges toward ~ 0.05, which is congruent with the value predicted in eq. j4.22l >: sin 2 6 ~ 0.05. 
For high densities (x < 20), matter suppression is severe. No flavor oscillations happen and 
the off-diagonal overlap functions fi ell and I efl are close to zero. In free space for ~x > 60, the 
imaginary part I efi oscillates symmetrically around zero whereas the real part % eii is positive. 



4.3.1 Equivalence with wavefunction quantum mechanics 



The solution given in Fig. l4.5l is numerically equivalent to the solution obtained using the stan- 
dard wavefunction formalism (Giunti 2004, Wolfenstein 1978, Bethe 1986, Mikheyev & Smirnov 
1986): 



(4.23) 



-A 




-\{ 














-1 







when we identify the specific intensities with the probability densities 



K " l«AvJ 2 

Iy u « l«AvJ 2 , (4.24) 



and the macroscopic overlap functions with 

leu " h ~ KK) ■ ( 4 - 25 > 



Thus, our Boltzmann formalism is completely consistent with the existing description. 
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4.4 Conclusions 

In this chapter, we have demonstrated that our generalized Boltzmann equations for the neutrino 
specific intensities, eq. J3.34> . reproduce the expected behavior. We have considered three differ- 
ent settings, for which we numerically solved dimensionless versions of eq. (13.341 . 

• The vacuum beam solution was calculated in analytic form for a non-interacting neutrino 
ensemble in a box, eq. ( 14. 8> . 

• We have shown how the interplay between flavor oscillation and decoherent matter coupling 
can be analyzed simply be solving eq. J4.7i for neutrinos undergoing absorptions on nucle- 
ons. The graphical results for three different regimes are given in Figs. 14. 1114.21 and !4.3l 

• An important feature of our new approach is that the MSW effect for a streaming neutrino 
beam can be analyzed in a very simple fashion for observable quantities, see eq. (I4.2H . Nu- 
merical results for an exemplary density profile are depicted in Fig. 14.51 

• For realistic physical settings in a supernova environment, one must take into account the 
energy spectrum of the neutrinos, energy redistribution scattering, etc., and calculate the so- 
lutions to eq. J3.34l > in their full generality. This has not been performed in this work and it 
would also be interesting to understand better the effects of the self-interaction contributions 
on the right-hand-side of eq. ( 13.27b . 
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Summary of Part I 



We briefly highlight the main features of our formalism in flat spacetime. The analytic expres- 
sions and conclusions reached in Part 1 are published in Physical Review D (Strack & Burrows 
2005). 

5.1 Quantum effects 

• Neutrinos are mathematically represented by spinor fields. One consequence of their spin-^ 
nature is the Fermi-Dirac state occupation probability. We include this in the stimulated ab- 
sorption terms of our formalism in eq. J3.36i . Furthermore, we argue that it is straightforward 
to include blocking corrections for the sink and source terms originating in the flavor oscil- 
lations in eqs. J3.27b and J3.34i as well as in the corresponding collision integrals (Burrows, 
etal. 2004). 

• In the simultaneous presence of collision and oscillations, two time scales govern the evolu- 
tion of the ensemble. When the collision or interaction time is of the same order as the flavor 
oscillation periods, it is believed that the collisions decohere the oscillation cycles and reset 
their clock. Then the temporal evolution becomes non-trivial since both processes must be 
taken into account for. We have shown in chapter|4]and in particular the in series of graphs, 
Figs. 14. l| 14.21 and 14.31 that our formalism encapsulates this interplay of flavor oscillations 
and matter coupling in a straightforward way. 

• The possibility of oscillation into a sterile neutrino species is explored in section l4~2l We 
can observe indirect affection of the sterile neutrinos, since their oscillating flavor partner, 
the v' e s do interact. The plots are shown in Fig. 14.41 This is potentially important for kinetic 
calculations in the early universe (Dolgov 1990) and could also have an effect on the flavor 
composition of the neutrino spectra emerging from supernovae (Abazajian, et al. 2004). 



56 5 Summary of Part I 



• Abstractly speaking, the neutrino vector is resonantly rotated in flavor space through mat- 
ter interaction with electrons (MSW effect) at certain densities. This has been suggested to 
explain the solar neutrino puzzle (Bethe 1986). We have embedded this effect in a simple 
quasi-classical framework and have shown analytically that our formalism correctly displays 
this important quantum-physical prediction in subsection 14.3.11 The corresponding plot is 
depicted in Fig. 14.51 



5.2 Classical framework 

• The classical framework and the classical equations of radiation hydrodynamics (Mihalas 
1999) are commonly employed to simulate core-collapse supernovae. The neutrinos are be- 
lieved to play a crucial role in the explosion process as they carry away almost all the gravita- 
tional energy. Future detection of neutrino spectra from supernova explosions will unravel the 
supernova mystery a little more and buttress or falsify the theoretical models on the market. 
Therefore it is very important to consistently include neutrino oscillations into the propaga- 
tion from the supernova to the detector but also to include it in the kinetic description close 
to the core. 

• Our generalized set of kinetic equations, eqs. J3.27I and J3.34l >. will allow straightfor- 
ward extension of any numerical code that simulates core-collapse, to account for the 
purely quantum-mechanical effects of flavor oscillations and self-interactions. The gener- 
alized Boltzmann equations without neutrino self-interactions, eq. J3.34i . are very simple 
and strongly resemble the classical equations of transport theory. 

• Conceptually, we have reduced the time evolution of neutrino field operators over Fock space, 
eq. ( 12.91 . to the time evolution of real numbers by ensemble-averaging, eq. J2. 16b . 

• We have connected three approaches: We started with quantum field operators, eq. ( 12. 9t , 
made the connection the equations for the quasi-classical phase-space densities and the as- 
sociated specific intensities in section IT. 3. II and finally we showed that our approach is also 
consistent with the one-particle wavefunction quantum mechanics in subsection l4.3.1l 



Part II 



Curved Spacetime 
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Neutrino Phenomenology 



In general astrophysical contexts such as the early universe and in the vicinity of supernova cores, 
one needs to consider the effects of gravitational fields on the physics under consideration. Here 
we seek to understand the kinetics of oscillating neutrinos interacting with a background medium 
in the presence of a strong, non-negligible gravitational field. 

Two changes to the eq. J3.27I for flat spacetime need be implemented. Firstly, we must ana- 
lyze a gravitational contribution to the mixing Hamiltonian (Cardall & Fuller 1997), which has 
observational consequences only if other effects such as spin precession in the presence of a 
magnetic field take place (Piriz, Roy & Wudka 1996, Briiggen 1998). This gravitational contri- 
bution stems from the spin connection, which we need to consider when treating spinor fields in 
a curved spacetime (Brill & Wheeler 1957, Weinberg 1972). 

In addition to the modification of the mixing Hamiltonian we must write the left-hand-side 
of the Boltzmann equations, eq. J3.27i . in a manifestly covariant manner so as to account for 
curvature. To this end we follow (Mihalas 1999 & references therein, and Stewart 1971). This 
will be done in section fTTI 



6.1 Neutrino oscillations (i) 

We discuss how mixed neutrinos can be analyzed in a covariant manner following (Cardall & 
Fuller 1997). We set G — ft = c — 1 throughout this section. In a standard treatment, the neutrino 
flavor state is written (Kayser 1981) 




(6.1) 



where flavor (mass) indices are in Greek (Latin letters), and the transformation coefficients are 
given by 
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U 



sin 9 cos 
- cos 6 sin 



(6.2) 



The mass eigenstates are taken to be energy eigenstates with common energy E. The three- 
momenta of the mass eigenstates are then 



/ m i 

> f = Je 2 -m 2 *E-—, 

1 y j 2E 



(6.3) 



where rrij is the rest mass corresponding to the mass eigenstate \mj). To observe oscillation it is 
crucial that one assumes null trajectory (x = t) of the mixed neutrino yielding the expression 



v a (x)) = ^ U a jS^p(-i^xj\mj) . 



(6.4) 



If the mass eigenstates could be measured at different positions (or times), the interference pattern 
would be destroyed. How can eq. ( 16.41 be generalized to curved spacetimes and be cast in a 
manifestly covariant form? 

The transformation coefficients from eq. ( 16. 21 (or equivalently the vacuum mixing angle 6 
between the two flavor species of interest), the phase and the mass eigenstates are frame inde- 
pendent quantities. The trajectory of a free neutrino is null geodesic, which leads to 



(6.5) 



In this expression, P 1 is the four-momentum operator that generates spacetime translation of the 
mass eigenstates. The quantity is the (null) tangent vector to the neutrino's world line, 



xf(A) 



xi(A) 
x 2 (A) 
x 3 (A) 



(6.6) 



and A is an affine parameter of the world line. In the following, we prove that eq. (16. 5> is indeed 
the proper expression for mixed neutrinos in a curved spacetime. 

Therefore, we consider motion along the x\ axis. Then the four-momentum operator reads 



P" = (£,P,0,0) , 
and with the Lorentzian signature metric 

ij mv = diag[-l, 1, 1, 1] 



(6.7) 



(6.8) 



we have for the argument of the exponential in eq. d6.5t 
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^dt dx 
dA 1 dA 



dA 



In this expression the X\ -component of the four-momentum operator is 

M 2 

Pl = Vl ^,E--. 



(6.9) 



(6.10) 



M is the mass operator which we replace by its eigenvalue and write for the corresponding phase 

CD : : 



, dt/dA 
dx/dA 



COj = -i f 

J™ m 2 
-r—dx 
* 2£ 



m 2 ^ 
E - 

2E 



^dA 
dA 



— —i — - (x — Xn) . 

2E 

Here we have used the null trajectory condition 

dt/dA _ p' null 
dx l dA Pitt 



= 1 



(6.11) 



(6.12) 



The important result is that eq. d6.5i is the covariant generalization of eq. J6.4i to spaces with 
curvature. 



6.1.1 Matter-induced effective mass for the v e 

We consider Dirac neutrinos and the equation of motion is the Dirac equation. We can include 
matter effects into this equation. For this purpose the Dirac equation in the flavor basis can be 
cast in the form, 

[ir(d fi + A fl P L ) + M]ifr = 0. (6.13) 

Here iff is a column spinor with four components for each flavor, y are defined by eq. fl2.3i and 
our choice shall be the chiral basis (Peskin & Schroeder 1995): 





(6.14) 



where cr' are the Pauli-matrices. Vl is the chirality projection on the left-handed components 
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l - r 5 

Pl = —^~, (6-15) 

where 

i% 

Here one makes the assumption that only the left-handed components of the flavor spinor ijj 
interact, whereas the right-handed components are taken to be sterile. The mass matrix in flavor 
basis M satisfies, 



mV 



M-t'ln 1 ^^, (6-17) 



where U has been specified in eq. ( 16. 2> and involves the vacuum mixing angle 0. 

A 11 is the flavor-basis effective potential matrix for an interaction with the electron back- 
ground: 

A" = -y/2G F ^ty . (6.18) 

In this expression N% is the number current of the electron fluid 

N% = n e u?, (6.19) 

where n e is the electron density in the fluid rest frame, and u 1 * is the fluid's four-velocity. This is 
the covariant extension of matrix given in eq. ( 13.301 . 

By iteration of the Dirac equation eq. (16. 1 3i in Fourier space one has for the mass shell 
relation the expression, 

(J* + A"P L ) (Pf, + AJP L ) = -M . (6.20) 
Keeping terms first order in Gf and employing eq. ( 16.71 we find 



2E 

with 



(6.21) 



V = -2^G,(^S) ■ (6-22) 

The three-momentum operator eq. I6.2H now includes effective mass contributions from back- 
ground matter, and can be used in eq. ( 16. 5> . 

The most important physical consequence of this effective mass is well known (MSW effect) 
and we have performed explicit calculations in this thesis in section l431 
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6.1.2 Spinor fields 

Spinor fields and their transformation properties are investigated with help of the vielbein for- 
malism which separates local Lorentz transformations and general coordinate transformations. 
This is convenient because for spinor fields one can then apply, locally, the spinor representation 
of the Lorentz group for which the well known y's of the Dirac algebra in flat spacetime eq. 
fl2.3i can be employed. In Appendix 151 we sketch how to construct locally inertial coordinate 
systems and "hide" the curvature in the components of the vielbeins so as to apply Lorentz trans- 
formations on spinors in curved spacetime (Carroll 2004, Weinberg 1972, and Brill & Wheeler 
1957). 

The Dirac equation in curved spacetime involves gravitational effects on the spin through the 
spin connection r^. 

(rfe a "[d ll + r il ]+M)tfr=0, (6.23) 

where Greek indices refer to the general coordinates, while the Latin indices a, b, c, d refer to the 
locally inertial noncoordinate basis as described in Appendix |Bj Using the explicit expression 
for r» (Weinberg 1972): 

r^ l -[^,y c \e b v e cv ,, , (6.24) 

where the semi-colon denotes the covariant derivative, one can incorporate the gravitational ef- 
fects on spin into the momentum operator P 1 * of eq. d6.5i . The Dirac equation J6.23t features the 
Dirac matrix product 

/ [/, f] = 27f V - 2/7°'/ - 2ie dahc y,y d , (6.25) 

where the right-hand-side can be calculated by using the identities for Dirac matrices as for 
instance given in (Peskin & Schroeder 1995). Here is the totally antisymmetric tensor with 
f oi23_j After some manipulation we can group the terms arising from the spin connection as 
follows: 

Wr M = W((S^(-^)) , (6.26) 

with 

®M = 5* Ube W 1 ydeb v e crv . (6.27) 

Note that in the reference (Cardall & Fuller 1997), a different convention of the y's is employed 
and the corresponding expression for ©^ is different (I could not see how the authors derived their 
expression). It is now apparent that, by adding a term proportional to the identity matrix (this has 
no observable consequence), we can choose either chirality to be prone to the gravitational effects 
and write 
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ryWr„ = r/e/ (ffi„P L ) , (6.28) 

where 

i-r 5 

= — ^- (6.29) 

is the projector on the left-handed states. 

An important point is that the gravitational contribution & is proportional to the identity ma- 
trix in flavor space, and diagonal in spin space. It cannot induce spin flips on its own. Therefore it 
will not have any observable effects unless there are other off-diagonal terms in spin space, e.g.: 
from the interaction of a neutrino magnetic moment with a magnetic field (Piriz, Roy & Wudka 
1996). In reality, however, the spin-gravity coupling needs a large anomalous magnetic moment 
of the neutrino to yield observationally relevant effects. 

At present the Standard Model-like calculation for massive Dirac neutrinos predicts (Shrock 
& Fujikawa 1977, Mohapatra & Pal 2003) 

3eGp 

m v , (6.30) 



8 V2tt 

where m v is the mass of the neutrino. In numbers, we have 

^3.1x10-'Vb(^) , (6.31) 

fiB being the Bohr magneton, ej (2m e ). This is rather puny. Even though the magnetic fields near 
protoneutron stars or supernova cores can be as large as 10 12 Gauss, the spin-flip length scale 
will be nonlocal, i. e. spread out over large spatial distances when compared with the scales of 
other processes as collision with the ambient matter. For theoretical reasons it is important that 
the gravitational contribution is added to the mixing Hamiltonian eq. fl3.3i even though with the 
present range of magnetic fields and anomalous magnetic moments this will have no physical 
effect. 

To find effects (Lambiase 2004, Lambiase, Papini, Punzi & Scarpetta 2005), the neutrino 
magnetic moment must be estimated to be magnitudes larger than that following from eq. J6.3H 
with neutrino masses in the eV range. 
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6.2 Neutrino oscillations (ti) 

We can now treat the effective mass contribution of eq. J6.18t and the spin connection of eq. 
J6.27> in the same fashion. Then the flavor state in curved space time of eq. d6.5i can be written 

as 

\v a (M = J] U aj exp |i ^ pvp^mj) . (6.32) 

In this expression the argument of the exponential now encompasses spin and effective mass 
effects and we have 

P M P M =-l^j-+P M [A tl + <6 lt ) Pl) ■ (6.33) 

We are in position to explicitly calculate the vielbeins eq. j6.24t for interesting geometries and 
follow the temporal evolution of the flavor composition of the mixed neutrino spinor in curved 
spacetime through solving the following modified version of eq. ( I6.23t : 

(iy a e/ + (A„ + ® M ) P L ] + M) ifr = . (6.34) 

More equations of this type are derived in the references Cardall & Fuller 1997, and Piriz, Roy 
& Wudka 1996. 



6.3 Conclusions 

• In this chapter, we have discussed how to generalize the phase of a mixed neutrino to curved 
spaces. The result is eq. J6.33i . 

• We have shown the dynamical equation j6.34t for a neutrino spinor in curved spacetime. 
We argued that the non-vanishing contribution due to spin-gravity coupling can be written in 
terms of the vielbein fields, eq. J6.27b . This contribution is diagonal in spin and flavor space. 

• The terms in the operator of eq. J6.34I . notably the gravitational contribution, (5^, will be 
used to write the mixing Hamiltonian, eq. il.161 . so as to account for spin-gravity coupling. 
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General-relativistic Boltzmann Formalism 



The purpose of section lTTl is to formulate the classical Boltzmann equation 
3/(r,p,f) . 3/(r,p,t) . df(r,p,t) 

^ 1" r ' ^ 1" P ' 7. 

dt dr F dp 

in a manifestly covariant manner and generalize the Liouville operator of the left-hand-side to 
curved spaces. In section 17121 we decompose the phase-space densities into their flavor degrees 
of freedom and into their spin degrees of freedom. We show the associated mixing Hamiltonian, 
eq. f7l6l . 

The highlight of this chapter are the new kinetic equations for oscillating neutrinos in curved 
spaces, eq. il. 241 . 



Sf 
St 



(7.1) 



7.1 Classical Boltzmann equation 

The phase-space for a system of relativistic particles is simply the tangent bundle over the space- 
time manifold, and the corresponding distribution function or phase-space density can naturally 
be described as a scalar field on this bundle (Stewart 1971). At every point x of the semi- 
riemannian manifold of our spacetime, the tangent space, T x , becomes the relativistic generaliza- 
tion of momentum space. Physically realizable four-momenta are timelike or null. Free particles 
have a geodesic trajectory and between collisions, the four-momenta of the particles, 

1A^> (7 ' 2) 
where A is an affine parameter, satisfy the geodesic equation: 

% + r\ pP v f = . (7.3) 

uA 
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= + (7-4) 



(7.5) 

coll 



In this expression r M Y „ is the symmetric, metric-compatible and unique Christoffel connection 
that is linear in the derivatives of the metric: 

2 8 \dxP dx v dx a 

where g MV is the twice-covariant metric tensor field which multi-linearly makes two tangent vec- 
tors to real numbers. 

The manifestly covariant relativistic Boltzmann equation for the invariant phase-space den- 
sity / (x^, p M ) is then given by: 

The term on the right-hand-side is the collision term (Mihalas 1999). 

There is a subtlety involved in this way of writing the Boltzmann equation. We have tacitly as- 
sumed that the invariant phase-space density / (xf , p**) is defined for all possible four- momenta. 
However only the four-momenta that lie on the null-cone make sense for the description of al- 
most massless neutrinos (as in chapter[6]we assume null geodesies for the neutrinos). Therefore, 
we must ensure that remains on the null cone as the neutrino propagates. This can be achieved 
by employing the vielbein formalism (AppendixlBl because in the local vielbein frame we can fix 
the four-momentum to be null once and for all, since in the vielbein frame we have the Minkowski 
metric. All the spatial-dependence involving the curvature is "out-sourced" to the complicated 
structure of the vielbeins. We regard the phase-space density / as a function of (xf,p a ), where 
the components of the spatial dependence are specified with respect to the global coordinate 
system. The momenta p", however, are given in the local vielbein frame. Then one has, in a kind 
of mixed representation similar to eq. (17.51 . 



,dl + dp^dl 
a dxf dA dp b 



Sf 
SA 



(7.6) 

coll 



fin 

We need to calculate Recalling that the neutrino trajectories are geodesies fulfilling eq. 
fl7.3i in the global coordinate system, there is 

Jo 



dp b ^ dp 

dpi 1 def c c 



dA ' 11 c dA 



i-' 



dA dA P 



= e b ll (-r» po .p°p'r-e»c, P P e p c ) 



-e\e£e» Cfi p a p c 

-^\ c p a p c . (7.7) 
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In this formula, the comma denotes a partial derivative and a semi-colon stands for the covariant 
derivative. We can define the so-called Ricci rotation coefficient 



9f 



and the Pfaffian derivative (this is the correct partial derivative in the vielbein frame): 

d a - e/dfi . 

In terms of these quantities the Boltzmann equation eq. J7.5I > becomes 



/ = 



Sf 

SA 



coll 



(7.8) 



(7.9) 



(7.10) 



This is the general form of the left-hand-side of the Boltzmann equation (sometimes referred to 
as the Liouville operator) in curved spacetime. In non-negligible gravitational fields this Liou- 
ville operator eq. J7. 101 replaces the flat space left-hand-side of the Boltzmann equation as for 
example shown in eq. (17. 1> . 



7.2 4 x 4 Mixing Hamiltonian 



We increase the dimensionality of the mixing Hamiltonian employed in flat spacetime, eq. 0.3I >. 
to further decompose the phase-space densities and the associated intensities not only into their 
flavor degrees of freedom, but also into their spin degrees of freedom. The matrix elements of 
the Wigner phase-space density are then written as a 4 X 4 matrix: 



5 = <";lpl"/> 

I A A Al Al\ 

Jv e Je[i Jv e Jefi 

A A Al Al 

J pe JVp J[ie JVp 

rlT rIT fl A 

Jv e Jefi Jv e Jefi 

rlT rIT fl A 

\Jue Jv u J fie Jv u l 



(7.11) 



In this matrix i and j run over flavor degrees of freedom and the spin states, f denotes left-handed 
states and J, means right-handed states. In addition to the previously employed macroscopic over- 
lap for the differently flavored neutrinos in the ensemble, eq. J2.20I I. we must also consider mixed 
spin states denoted by the superscripts fl and it- Again this is only relevant if the magnetic mo- 
ment of the neutrino is comparatively large and the magnetic fields are very strong (Piriz, Roy & 
Wudka 1996). The interaction of the neutrino spinor with the electromagnetic field stems from 
terms of the form 



(7.12) 
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where i/r is the neutrino spinor, as for example defined in eq. ( I6.23> . and yU the magnetic moment 
of the neutrino, eq. ( 16. 3U . For S ah we have the expression eq. ( IB.9I from the AppendixlBl 

S ah = l -[y\y h \ . (7.13) 
This electromagnetic interaction term can be added to the Dirac equation J6.34i . yielding 

{i-fef [d M + (A M + «„) P L ] + ftS ab F ab +M)tfr=0. (7.14) 

In the mixing Hamiltonian, the surviving electromagnetic terms are [iB with B being the com- 
ponent of the magnetic field perpendicular to the neutrino trajectory. Similar to as performed in 
the papers by Piriz, Roy & Wudka 1996, and Briiggen 1998 (these references decomposed the 
neutrino wavefunction into their spin and flavor states), we generalize the mixing Hamiltonian 
of eq. J3.3I > as follows: 

2B(^, if) = 2B vac + + atU/A *") ~ 2Bw(^ , *?) + 2Bgrav(*") • (7.15) 

This Hamiltonian includes the weak, electromagnetic and the gravitational interactions. Adding 
up all the contributions, except for the self-interactions, the most general form of the mixing 
Hamiltonian is: 



<V •„ + <V ' X +P ll & t Is cos 20 ^ sin 20 ^B ^B ' 

sin 29 <V„ c +P II W/ s +d£ cos2 g ^ e B n v B 

fi Vr B fx eli B cos 26 4g sin 20 

/V B ^B % sin 29 4*L cos 29, 

(7.16) 



*y cc is the charged-current contribution leading to the effective mass of the v e in matter multi- 
ply discussed in this work (the covariant formulation is done in subsection l6.1.1l and the standard 
way is treated via the matter mixing Hamiltonian in eq. (13. 8» : 

'Vcc =2y/2G f n e . (7.17) 

<V nc is the neutral-cuiTent contribution (Mohapatra & Pal 2003): 

<V„ c =-2^ f n n , (7.18) 

where n n is the neutron number density and n e is the electron density. We emphasize that the 
neutral-current contribution is the same for all flavors of neutrinos (which is why it was dropped 
in eq. J3.8i . whereas the charged-current contribution affects the v e only (see Fig. 13. IV ). For 
neutral-current scattering on the neutrons in the ambient matter, we have the diagrams depicted 
in Fig. O 
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Fig. 7.1. Neutral-current scattering with neutrons for v c 's and v„'s. 



7.2.1 Generalized resonances 

By virtue of the mixing Hamiltonian eq. ( I7.16> we realize that resonances occur whenever the 
diagonal elements cross. We read off the resonance conditions for spin-flips and MSW conver- 
sions. 



For fl^fl 



for fl^fi: 



for fl^fi- 



Am 2 

^Vcc cos6» = 

2s 



'V cc + <V nc +P fl ®'/s = , 



Am 

<V CC + <Vnc + P^ls — — cos 26 — 
2s 



Am 

Vcc + *V„ C + PMle + — — cos 26* = 
2s 



T I 

and lastly for the transition /„' <-> f; : 



<V nc + P fl &'/ £ = 



(7.19) 



(7.20) 



(7.21) 



(7.22) 



(7.23) 



In the active-active MSW transition, eq. J7.19> . the matter potential is due to the charged-current 
interactions only and proportional to the electron number density (see eq. fl7. 171 . To obtain ac- 
curate spectra in general conditions all these resonances must be taken into account. The grav- 
itational contribution generates a shift of the resonance positions. The adiabaticity can also be 
altered (Bruggen 1998). 
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7 General-relativistic Boltzmann Formalism 



7.3 Kinetic equations 

With the aid of eq. (17. 10i . the associated generalization to curved spacetime including the spin- 
states of the the Heisenberg-Boltzmann equation, eq. P. Il l, reads: 



This equation is the main result of Part II of this thesis. Note that g and 2C are now 4x4 ma- 
trices. An analogous equation can be written for antineutrinos. Particles and antiparticles couple 
identically to the gravitational field, so no changes need to be made for the gravitational modi- 
fications; neither for the Liouville operator in curved space, eq. ( I7.10t . nor for the gravitational 
contribution © in eq. J6.27t to the mixing Hamiltonian 9B, eq. d7 . 1 61 . 

If we were to componentize them like done above in eq. P.27l >. we would have to deal with 
sixteen coupled equations for the neutrinos without the self-interactions. The issue of making 
the equations real-valued can naturally be solved by defining the real and imaginary parts of the 
off-diagonal elements of g similarly to eq. J3.22l >. Including the self-interactions, however, one 
would have to solve 32 nonlinearly coupled equations. Clearly, for real applications one would 
have to break up the neutrino trajectory to care for the most dominant effect one at a time. 

Eq. J7.24I can also be used if one wishes to generalize eq. J3.27i to curved space without 
taking to account the coupling of the still speculative anomalous magnetic moment to a very 
large magnetic field. Then g and 2C reduce to the familiar 2x2 matrices given in eqs. ( 12. 161 and 
fl3.3t . In this case the only effect of the curvature is the standard generalization of the Liouville 
operator to curved space as performed in eq. (17.101 1. 



7.4 Conclusions 

• The basic fundament of eq. 17 . 241 is the classical, general-relativistic framework with which 
we generalized the Liouville operator to curved spaces in eq. J7 . 1 Oi . 

• The quantum physics is incorporated by means of the general mixing Hamiltonian, eq. (17. 16t . 
We have included spin degrees of freedom. 

• The "master equation" of Part II, eq. J7.24> describes, in principle, the evolution of the quasi- 
classical density in phase-space for massive neutrinos with collisions including all known 
effects - from the quantum-mechanical viewpoint as well as from the gravitational perspec- 
tive. 




m, g]+c. 



(7.24) 
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Summary of Part II 



We briefly highlight the main features of our formalism in curved spacetime and discuss physical 
applications for which this formalism might be helpful. 

8.1 Spin-gravity coupling 

Invariance under general coordinate transformations for spinors naturally leads to the vielbein 
fields introduced in AppendixlBl In short, these considerations lead to an additional contribution 
for the differential operator of the Dirac equation J6.23l > in terms of the spin connection. How- 
ever, terms arising from the spin connection are proportional to the identity matrix in flavor space 
and do not lead to additional effects concerning oscillation between flavor species. Moreover, the 
spin connection terms are also diagonal in spin space and act on one spin orientation of the neu- 
trino only. 

The most general mixing Hamiltonian, eq. J7.16> . has potentially other off-diagonal contri- 
bution stemming from spin-flips from coupling of the neutrino spinor the an external magnetic 
field. It is for this case only that the gravitational resonances in subsection l7.2.1l are important. 
It has been analyzed elsewhere that, assuming large anomalous magnetic moments, these effects 
can be potentially large (Piriz, Roy & Wudka 1996). 

8.2 Possible applications 

Combining the oscillation physics in curved spacetime from section [6] with the kinetic general- 
izations from section f7~Tl vields the master equation of Part II, eqs. j7.24> . These equations satisfy 
all quantum-mechanical and gravitational considerations for the kinetics of oscillating neutrinos 
in curved spacetime. 
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Potential applications include active galactic nuclei (Piriz, Roy & Wudka 1996) and the neu- 
trino fluxes in supernova cores. In practice, however, one would break up the master equations, 
eq. J7.24I I. and solve them isolating the dominant effect. For supernova neutrinos, for instance, 
we do not need to care about collisions at nuclear densities because observable flavor oscillations 
are matter-suppressed. There, we must focus on the collision physics with the ambient matter. 
In the streaming region above the neutrino sphere, self-interactions become important and they 
must be focussed on. Spin-flips only deserve consideration for very large magnetic fields and so 
forth. 

To close this summary of Part II, we note that it has been proposed that the nonlinear evo- 
lution in a neutrino background leads to a speed-up of flavor equilibration processes that could 
potentially overshadow the gravitational and all other effects (Saywer 2004 and 2005). We leave 
this question open for future work. 



Part III 



Conclusions 
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Conclusions 



In this thesis, we have derived a generalized set of Boltzmann equations for real-valued phase- 
space densities of oscillating neutrinos interacting with a background medium. The off-diagonal 
functions of the Wigner phase-space density representing macroscopic overlap are explicitly in- 
cluded and serve to couple the flavor states to reflect neutrino oscillation physics. Conceptually, 
we have reduced the time evolution of creation and annihilation operators to that of real val- 
ued phase-space densities without losing quantum-physical accuracy. Important quantum effects 
such as matter-enhanced resonant flavor conversion and quantum decoherence through matter 
coupling are correctly incorporated. The generalized Boltzmann equations are simple and very 
similar to the equations of classical transport theory. Neutrino oscillations are incorporated by 
new sink and source terms that indirectly couple the expanded set of equations. Using this for- 
malism, codes that now solve the standard Boltzmann equations for the classical neutrino phase- 
space density (/„.), or which address its angular and/or energy moments, can straightforwardly 
be reconfigured by the simple addition of source terms and similar transport equations for over- 
lap densities that have the same units as / v ., to incorporate neutrino oscillations in a quantum- 
physically consistent fashion. In curved spacetime, oscillation phenomenology is altered. These 
effects play no observable role except for the usual redshift and the contribution to the mixing 
Hamiltonian from the spin connection can be isolated in the diagonal of the left-handed spin 
sector of the effective mixing Hamiltonian. If, however, spin flips occur through coupling of an 
anomalous magnetic moment of the neutrino with a large magnetic field, one must take into 
account gravitational resonances similar to the MSW effect in flat spacetime. 



Part IV 



Appendices 



A 

Cross section v e + n — > e~ + p 



A convenient reference neutrino cross section is cr , given by 

1.705 x cm 



cr = 



4G2(m e c 2 ) 2 „ nc „ 1A _ 44 2 



n(ficy 



(A.l) 



where is the Fermi weak coupling constant 1.436 x 10 49 ergs cm 3 ). The total v e — n 
absorption cross section for the reaction v e + n — > e~ + p is then given by 



2 \ 



np 



1 - 



m e c 



1/2 



(A.2) 



where gA is the axial-vector coupling constant (~ -1.23), and A np - m n c 2 - m p c 2 = 1.29332 
MeV . In Fig. (IA. li we show the heuristic Feynman graph corresponding to this process. 




Fig. A.l. Super-allowed charged-current v e absorption on neutrons. 



B 

Vielbein (Tetrad) Formalism and Spinor Fields 



The essential quantities in the general theory of relativity such as the curvature or the energy- 
momentum are invariant under general coordinate transformations. This can be viewed as the 
construction principle of the general theory of relativity (Weinberg 1972). Luckily, these quan- 
tities behave like tensors not only under general coordinate transformation but they also behave 
as tensors under Lorentz transformations. This makes the transition from special relativity in fiat 
spacetime to general relativity in curved spacetimes straightforward. For spinor fields this is not 
true because the tensor representations of the group GL(4) of general linear 4x4 matrices be- 
have like tensors under the subgroup of Lorentz transformations, but there are no representation 
of GL(4) which behave like spinors under the Lorentz subgroup. To remedy this situation, one 
erects at every point x a noncoordinate basis % a (x) that is locally inertial. Then the metric can be 
written 

gfiv(x) = e fl a (x)e Y b (x)7] ab , (B.l) 

where for a Lorentzian spacetime rj a b represents the Minkowski metric, while in a space with 
positive-definite metric it would represent the Euclidean metric. The Latin indices indicate that 
the locally inertial coordinates are not related to any coordinate system. The vielbeins are given 
by 

Coordinate tensors of arbitrary rank can be made into a set of scalars by expressing their compo- 
nents with respect to the locally inertial basis, that is, 

T a b = e%e b v T\ . (B.3) 



The point is that now, just as in global Minkowski space, the change of the locally inertial basis 
in the vielbein frame corresponds to a (now spatially dependent) Lorentz transformation. The 
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B Vielbein (Tetrad) Formalism and Spinor Fields 



two transformations, local Lorentz transformations and general coordinate transformations, can 
be executed at the same time resulting in the mixed transformation law (Carroll 2004): 

r'^ v = A"' fl |^A^— r^ v . (B.4) 
ox^ ox y 

In general an arbitrary field in the vielbein frame if/„(x) transforms as 

4r„{x) = J] [D (A{x))] nm <!s m (x) , (B.5) 

m 

where D(A) is a matrix representation of the Lorentz group tailored to the nature of the field. 
Heuristically, we have stored away the curvature in the vielbeins so as to now being able, for a 
spinor field, to apply the spinor representation of the Lorentz group. 

We now need a derivative operator, Xs a , that is a scalar with respect to the general coordinate 
transformations and a Lorentz vector with respect to local position-dependent Lorentz transfor- 
mations with matrix representation A a b (x). This ensures, that any action which depends on the 
various fields of interest, including spinor fields, will then automatically be independent of the 
choice of locally inertial frames if it is invariant under ordinary constant Lorentz transformations. 
The derivative operator must transform homogeneously: 

D fl (A(x) -> A a b (x)D (A(x)) D b iff(x) . (B.6) 

It can be shown that the expression following expression is the correct derivative operator: 

= e/ + r,) , (B.7) 

where the spin connection r M is given by (Weinberg 1972): 

r p = \s bc e b v e cvv , (B.8) 

where the semicolon denotes the covariant derivative, and S hc is the appropriate representation 
of the Lorentz group. For spinor fields, employing the chiral choice of y's, eq. ( I6.14t . one has 

s 6c = ^[rV], (B.9) 

where square brackets denote the usual matrix commutators. 

To summarize: The effects of gravitation on any physical system can be taken into account 
by writing down the matter action or the field equations that hold in special relativity, and then 
replacing all derivatives d/dx a with the derivative operator £>„ of eq. iB.6\ . For the case of consis- 
tently describing spinor fields in curved space one must resort to this line of thought, since only 
in the vielbein frame can one employ the spinor representations of the Lorentz group coherently. 



References 



Abazajian, K., Bell, N., Fuller, G. M., and Wong, Y. Y. Y: astro-ph/0410175 (2004). 
Abazajian, K., Fuller, G. M., and Patel, M.: Phys. Rev. D 64, 023501 (2001). 

Akhiezer, A. I., and Peletminskii, S. V.: Methods of Statistical Physics (Pergamon Press, Volume 104, 1981). 
Akhmedov, E., Barroso, A., and Keraenen, P.: Phys. Lett. B 517, 355 (2001). 
Bell, N. F, Rawlinson, A. A., and Sawyer, R. F: Phys. Lett. B 573 86-93 (2003). 
Bethe, H.: Phys. Rev. Lett. 56 1305 (1986). 

Bilenky, S. M„ and Petcov, S. T.: Rev. Mod. Phys. 59, 671 (1987). 
Bilenky, S. M., and B. Pontecorvo, Phys. Lett. B 61, 248 (1976). 
Boyanovsky, D., and Ho, C. M.: Phys. Rev. D 69, 125012 (2004). 
Brill, D. R, and Wheeler, J. A.: Rev. Mod. Phys. 29, 3 (1957). 
Bruggen, M.: Phys. Rev. D 58, 083002 (1998). 
Burrows, A., Gandhi, R.: Phys. Lett. B 243, 149 (1995). 

Burrows, A., Young, T., Pinto, P., Eastman, R., and Thompson T. A.: Astrophys. J. 539, 865 (2000). 

Burrows, A., Reddy, S., and Thompson, T. A.: astro-ph/0404432 (2004). 

Cardall, C. Y, and Fuller, G. M.: Phys. Rev. D 55, 7960 (1997). 

Carroll, S. M.: Spacetime and Geometry, 1st edn (Addison Wesley 2004). 



86 References 



de Groot, S., van Leeuwen, W., and van Weert, C. G.: Relativistic Kinetic Theory (North Holland, 1980). 

Dolgov, A., Barbieri, R.: Nucl. Phys. B 349, 743-753 (1990). 

Friedland, A., and Lunardini, C: JHEP10 043 (2003). 

Friedland, A., and Lunardini, C: Phys. Rev. D 68, 013007 (2003). 

Fukuda, S., etai: Phys. Lett. B 539, 179 (2002). 

Giunti, C: Neutrino Mixing and Oscillations, www.nu.to.infn.it/ (2004). 
Lambiase, G.: astro-ph/04 12408 (2004). 

Lambiase, G., Papini, G., Punzi, R. and Scarpetta, G.: gr-qc/0503027 (2005). 

Kayser, B.: Phys. Rev. D 24, 1 10 (1981). 

Landau, L. D.: Sov. Phys. JETP 5, 101 (1957). 

McKellar, B. H. J., Thomson, M. J.: Phys. Rev. D 51, 4 (1992). 

Mihalas, D., and Weibel-Mihalas, B.: Foundations of Radiation Hydrodynamics (Dover Publications, 1999). 

Mikheev, S., and Smirnov, A. Y.: Sov. Phys. JETP 64, 4 (1986). 

Moertsell, E., Bergstroem, L., and Goobar, A.: Phys. Rev. D 37, 1237 (2002). 

Mohapatra R. N, and Pal, P. B.: Massive Neutrinos in Phys. and Astrophys. (World Scientific, 2003). 

Pantaleone, J., and Kuo T. K.: Rev. of Mod. Phys. 61, 937 (1989). 

Pantaleone, J.: Phys. Lett. B 287, 128 (1992) A. 

Pantaleone, J.: Phys. Rev. D 46, 510 (1992) B. 

Pantaleone, J.: Phys. Lett. B 342, 250-256 (1995). 

Parke, S. J.: Phys. Rev. Lett. 57, 1275 (1986). 

Peskin, M. E., Schroeder D. V.: An Introduction to Quantum Field Theory (Westview Press, 1995). 
Piriz, D., Roy, M. and Wudka J.: Phys. Rev. D 54, 1587 (1996). 
Qian, Y.-Z., and Fuller, G: Phys. Rev. D. 20 (1995). 



References 



Raffelt, G.: Stars as Laboratories for Fundamental Physics (The University of Chicago Press, 1996). 

Raffelt, G., and Sigl, G.: Nucl. Phys. B 406, 423 (1993). 

Raffelt, G., and Sigl, G: Astropart. Phys. 1, 165 (1993). 

Raffelt, G., and Stodolsky, L.: Phys. Rev. D. 37, 1237 (1988). 

Raffelt, G, Sigl, G., and Stodolsky, L.: Phys. Rev. D 45, 1782 (1992). 

Raffelt, G, Sigl, G., and Stodolsky, L.: Phys. Rev. Lett. 70, 16 (1993). 

Robinson, D. C: Phys. Rev. Lett. 34, 901 (1975). 

Rudszky, M. A.: Astrophys. Space Sci. 165, 65 (1990). 

Sawyer, R. R: Phys. Rev. D 42, 3908 (1990). 

Sawyer, R. R: hep-ph/0408265 (2004). 

Sawyer, R. R: hep-ph/0503013 (2005). 

Shrock, R. E., and Fujikawa, K.: Phys. Rev. D 16, 1444 (1977). 

Stewart, J. M.: Non-Equilibrium Relativistic Kinetic Theory (Springer- Verlag, 1971). 

Strack, P., and Burrows, A.: Phys. Rev. D 71, 093004 (2005). 

Takahashi, K., Sato, K., Burrows, A., and Thompson, T. A.: Phys. Rev. D. 68, 113009 (2003). 
Thompson, T. A., Burrows, A., and Horvath, J. E.: Phys. Rev. C. 62, 035802 (2000). 
Weinberg, S. W.: Gravitation and Cosmology (Wiley, New York, 1972). 
Wigner, E.: Phys. Rev. 40, 749 (1932). 

Wigner, E., Hillary, M., O'Connell, R. R, and Scully, M. O.: Phys. Rep. 106, 121 (1984). 
Wolfenstein, L.: Phys. Rev. D 17 (1978). 
Wolfenstein, L.: Phys. Rev. D 20, 2634 (1979). 



